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I.  INTRODUCTION 


It  has  been  recognized  lor  some  time  thai  one  of  the  major  shortcomings  of  current 
combustion  instabilit>  prediction  technology  is  that  it  is  computationally  removed  from 
conventional  interior  ballistic  analysis,  t  raditionally,  interior  ballistics  data  have  been 
extracted  from  an  analy  sis  and  used  to  formulate,  by  hand,  the  problem  to  be  run  with  one  of 
the  existing  combustion  instability  prediction  finite  clement  codes.  In  an  effort  to  initiate  a 
remedy  for  this  inefficient  operation  and  thereby  develop  the  basis  fora  more  elaborate  and 
meaningful  design  procedure  for  solid  propellant  rocket  motors,  the  computer  code 
CiRNMSH  (Grain  Mesh  I  was  developed.1  This  code  provided  an  automated  procedure  for 
generating  the  input  to  a  finite  element  mesh  generator  from  the  output  of  an  interior  ballistics 
code.  The  code  CiRNMSH  was  designed  to  couple  the  output  of  the  5(*4  Interior  Ballistics 
Computer  Program,  developed  by  Aerojet-General  Corporation,  with  the  computer  code 
FI  FSH.t  (Fluid  Mesh,  Generation.  2  Dimensions)  which  was  developed  as  a  companion  finite 
clement  mesh  generation  program  to  the  combustion  instability  prediction  code  FLAP.) 
(Fluid ^Analysis ^Program. ^Dimensions).*  The  development  of  GRNMSH  was  considered  to 
be  of  a  “proof  of  concept"  nature  and  stopped  short  of  pros  iding  the  complete  coupling  of  the 
two  codes  for  the  entire  period  of  performance,  i.e..  for  the  duration  of  propellant  burning. 
The  GRNMSH  code  generates  the  input  to  FI  I- SHI  only  at  the  instant  of  ignition.  This  report 
is  concerned  with  the  extension  of  the  GRNMSH  code  to  enable  the  same  automatic  transfer 
of  data  between  the  two  codes  at  any  instant  in  time,  from  propellant  ignition  to  burnout.  The 
extension  of  GRNMSH  has  resulted  in  the  development  of  the  code  BRNMSH  (Burn  Mesh) 
which  generates  the  input  to  FI  FSH.t  and  thus  initiates  the  generation  of  a  coupled  fluid-solid 
three-dimensional  finite  element  mesh  for  any  specified  depth  of  burn. 

The  development  has  concentrated  on  the  star  design,  which  is  depicted  in  Figures  I  and  2. 
Similar  designs  can  easily  be  handled  in  the  same  manner.  A  brief  discussion  of  the  previously 
developed  codes,  along  with  a  description  of  BRNMSH  follows. 


1  R  M  Hack'll.  4  CourM  Interior  Ballmicr- f'iniir  flemrni  i  omhusti,*  Insubilm  Analysis  Fr,\r,lurr  OS  Army  Missile 
Research  and  IVselopment  Command.  Redstone  Arsenal.  Alabama,  t  echnical  Report  1'- 78-72.  July  1*78 

2  X  I  W  hetstone,  I  R  three*  it.  and . I  S.  Hillheimer.  Hum,  Oram  IVogn  ami thrSb*  Interior  Hnlluiies  Computer  mm. 

St  M-UJ.  Interior  Ballistics  ITcpartment,  Applied  Mechanic*  Division.  Aerotet-Cieneral  Corporation.  Sacramento,  California, 
lunc  1^61. 

'  R  M  Haclett  and  R  S.  Juruf.  -I  three- Ihmemi,  mo/  finite  Clement  <  , s/e  for  Combust  ton  Insinbilm  PreJietnm 
Proceeding*  ol  the  I.Xth  JANNAF  Combustion  Meeting.  Nasal  Postgraduate  School.  Monterey.  California.  CPI  A  Publication 
281.  December  l*70 
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2.  THREE-DIMENSIONAL  FINITE  El  EM  ENT  MESH  GENERATION 


I  he  El.  A  P3  code  performs  u  linear  acouslo-modal  analysis  ol  the  irrotational  motions  of  an 
inviscid.  compressible  fluid  coupled  to  the  motion  ot  a  nearly  incompressible,  linearly 
viscoelastic  solid,  and  a  linear  potential  How  analysis  of  the  irrotational  motions  of an  inviscid, 
incompressible  fluid,  and  then  determines  the  effect  of  the  flow  field  and  of  combustion  on  the 
calculated  acoustic  oscillations.'  This  combustion  instability  analysis  is  performed  at  different 
points  in  time,  beginning  at  ignition,  or  time  zero.  The  output  of  FLAP3  is  modal  frequency 
and  an  evaluation  of  the  pressure  growth  decay  coefficient  for  each  mode  of  vibration 
analyzed.  A  positive  net  value  of  the  coefficient  indicates  a  growth  of  pressure  oscillations  and, 
therefore,  instability  while  a  negative  value  of  the  coefficient  is  an  indication  of  decaying 
oscillations,  or  stability. 

The  FI.  AP3  analysis  utilizes  the  finite  element  method  and  models  the  acoustic  cavity  and 
propellant  grain  as  an  assemblage  of  three-dimensional  quasi-hexahedral  (each  hexahedron  is 
a  combination  of  five  tetrahedral  elements  connected  at  the  corners.  In  order  to  be  realistically 
applicable,  any  three-dimensional  finite  element  code  requires  a  companion  mesh  generator;  a 
coupled  fluid-solid  analysis  code  requires  a  mesh  generator  of  more  than  the  customary 
complexity.  The  code  FLESH3  was  designed  to  be  used  in  conjunction  with  FLAP3.  A 
detailed  description  of  the  use  of  FLESH3  is  found  in  Hackett’s  User's  Manual  for  FLAP3 .* 
The  code  FLESH3  generates  the  numbered  nodal  points  and  their  coordinates  and  identifies 
each  node  as  to  whether  it  lies  in  the  acoustic  cavity,  on  the  cavity-grain  interface,  or  in  the 
grain;  generates  the  quasi-hexahedral  elements,  designated  by  the  eight  numbered  nodal 
points  defining  each  element,  and  identifies  each  element  as  to  whether  it  is  a  cavity  element,  a 
cavity  element  adjacent  to  the  cavity-grain  interface,  or  a  solid  propellant  element;  and 
generates  the  cavity-grain  interface  element  surfaces  (burning  surfaces)  and  identifies  each  as 
to  the  direction  of  the  surface  normal.  This  information,  along  with  a  few  additional  input 
parameters  relative  to  gas  and  propellant  properties  and  type  of  acoustic  mode(s)  to  be 
analyzed,  is  necessary  and  sufficient  to  activate  a  combustion  instability  analysis  by  FLAP3. 
given  adequate  computational  facilities. 

The  input  to  F1.ESH3  is  in  the  form  of  a  set  of  global  curves,  defining  cavity  and  grain 
boundaries;  a  designation  of  a  sequence  of  points,  as  to  whether  they  are  cavity,  interface  or 
grain  points;  a  designation  of  blocks  of  elements  to  be  generated,  by  the  nodal  indices  of  the 
blocks;  and  a  designation  of  the  number  and  location  of  longitudinal  cross-sections, 
establishing  the  number  of  layers  of  elements. 

4.  R.  M.  Hackelt.  I’srr's  Manual  for  FLAPS,  US  Army  Missile  Research  and  l>evelopment  Command,  Redstone  Arsenal. 
Alabama,  technical  Report  t'K-77-4,  July  1977. 
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The  object i\c  is  that  ol  minimizing  the  amount  of  time  required  to  prepare  the  input  for 
H  AP3  and  the  likelihood  ol  crior  in  the  preparation  ol  that  data.  Through  the  use  of 
BRNMSH.  which  vull  be  subsequently  described,  all  of  the  necessary  geometry  and 
identification  input  data  tor  I  I  AIM  can  be  generated  with  the  input  of  seven  geometry 
parameters,  tor  the  star  design,  plus  one  additional  combustion  parameter. 

3.  INTERIOR  BALLISTICS  CODE 

I  he  interior  ballistics  analysis  code  used  in  the  formulation  of  BRNMSH  was  developed  by 
Aerojet-General  Corporation.  I  he  BRNMSH  geometrical  input  parameters  arc.  therefore, 
those  presented  in  that  code,  but  similar  parameters  would  be  obtained  from  a  consideration 
of  a  comparable  interior  ballistics  code.  The  description  of  BRNMSH  will  be  related  to  the 
Aerojet  code,  but  the  procedure  followed  in  the  development  is  general. 

4.  DEVELOPMENT  OF  THE  BRNMSH  CODE 

A  listing  of  the  formulated  computer  code  BRNMSH  is  found  in  Appendix  A.  The  code  is 
based  upon  the  star  design  found  in  the  report  by  Whetstone  et  al.'  Due  to  the  dihedral 
symmetry  provision  in  FLAP3.  only  the  smallest  repeating  segment  need  be  analyzed.  This 
segment,  for  the  general  star  geometry,  is  shown  in  Figures  I  and  2;  with  the  independent 
geometry  parameters  identified.  The  only  difference  between  Figure  I  and  Figure  2  is  that  the 
angle  <t>  is  zero  in  Figure  2.  This  condition  necessitated  a  slightly  different  internal  provision  in 
BRNMSH.  Figures  3  and  4  define  the  different  burning  zones  and  the  regression  of  the 
interface  surface  with,  again,  the  difference  in  the  two  figures  being  the  value  of  the  angle  <f>. 
Figure  3  shows  a  typical  finite  element  mesh  for  burning  in  Zone  I .  The  cross-hatched  portion 
of  the  segment  represents  the  grain  elements.  The  mesh  is  seen  to  be  similar  to  that  which 
would  be  generated  for  the  case  of  zero  burning.  Figure  b  shows  a  typical  mesh  for  burning  in 
Zone  2.  It  can  be  noted  that  the  generated  mesh  for  Zone  2  burning  differs  distinctly  from  that 
for  Zone  I  burning.  Figure  7  shows  a  typical  mesh  for  burning  in  Zone  3.  As  the  burn  depth 
increases  it  can  be  seen  that  the  cavity  geometry  becomes  more  simply  defined  and  this  fact  is 
thus  reflected  in  the  generated  meshes  of  the  changing  cavity-grain  cross-section. 

A  working  knowledge  of  F1.ESH34  is  necessary  in  order  to  understand  the  intricacies  of 
BRNMSH.  and  to  be  able  to  modify  and  add  to  it.  but  not  in  order  to  use  it.  An  understanding 
of  the  use  of  BRNMSH  can  be  enhanced  through  the  consideration  of  an  example.  Table  I 
lists  the  twelve  cases  of  the  star  design  upon  which  the  formulation  of  BRNMSH  is  based. 
Design  No.  D-6c  is  chosen  for  the  example.  Three  different  burn  depths  are  selected.  They 
arc:  0.05  in.  (burning  in  Zone  1).  1.50  in.  (burning  in  Zone  2).  and  3.00  in.  (burning  in  Zone 
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Figure  1.  Cross-section  of  repeating  segment  of  star  design  showing 
design  parameters  N,  Rc,  W,  Ya,  a,  L  and<t>. 


Figure  2.  Cross-section  of  repeating  segment  of  star  design  showing 
design  parameters  N,  Rc,  W,  Y*.  a  and  L;  ♦=  0. 


Figure  3.  Cross-section  of  repeating  segment  of  star  design  showing  the 
four  progressive  burning  zones. 


Figure  4.  Cross-section  of  repeating  segment  of  star  design  showing  the 
four  progressive  burning  zones;  <t>  =  0. 
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Figure  5.  Finite  element  mesh  for  repeating  segment  of  star  design  with 
Zone  1  Burning;  <t>  =  0. 


Figure  6.  Finite  element  mesh  for  repeating  segment  of  star  design  with 
Zone  2  Burning;  ♦  =  0. 


Figure  7.  Finite  •foment  m**h  tor  repeating  segment  of  star  design  with 
Zone  3  Burning;  «t»  ■  0. 


'I  I  ach  case  icqunes  the  input  ol  two  data  cards  to  RRNMSH  Ihf  Inst  cant  contains  the 
seven  paiuinelcis  Y  K..  \N  .  </*.  \  ...  o.  and  l  tin  that  order).  1  he  second  card  contains  the 
hm  n  distance.  HIR\  I  he  output  lot  the  thiee  sepaiate cases,  which  isthe  input  to  I  I  l-SII.I. 
is  shown  in  Appendix  It 

It  can  he  noted  that  the  length  ol  the  grain,  and  provision  lot  live  equidistant  sections  arc 
conditions  which  aie  mternallv  set  in  the  program,  hut  could,  alternatively,  he*  input  as  an 
additional  set  ol  paiaineters  It  can  also  he  noted  that  the  value  of*.  the  necessaiv  oflset  lavlius 
ol  the  nodes  closest  to  the  axis  origin,  is  mternallv  set  at  0  1  in 

Iheie  aie  othei  modifications  which  could  he  made  to  HKNMSII  in  order  to  increase  its 
tlesihilnv,  lot  instance  that  of  developing  the  capability  ot  generating  a  mesh  lot  a  tape  ted 
longitudinal  profile.  Possible  modifications  such  as  this  are  obvious,  but  othei  more  subtle 
changes  could  also  be  made  which  would  increase  the  generality  of  the  code 

V  <.'ONU  UNIONS 

\  computer  code  has  been  developed  which  couples  an  interior  ballistics  analysis  code  w  ith  a 
combustion  instability  analysis  code  and  facilitates  the  efficient  lormulatton  of  the  three- 


dimensional  finite  element  model  for  instability  analysis  at  any  designated  point  in  time  during 
the  entire  performance.  Thus,  the  rocket  motor  designer  is  provided  with  the  potential  for  a 
more  inclusive  and,  therefore,  superior  design  approach. 
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TABLE  1.  STAR  DESIGN  PARAMETERS-SYMMETRY  7 

<N=7) 


STAR 

(Design  No.) 

ra 

w 

in. 

a 

deg 

IS 

V. 

In. 

H 

L 

in. 

D-5a 

10 

35 

80 

00 

■a 

0075 

D-5b 

10 

35 

85 

00 

BIS 

D-5c 

10 

35 

75 

□ 

0.136 

D  6a 

10 

35 

85 

0.0 

05 

0092 

4.1 

D-6b 

10 

35 

80 

00 

0  118 

D-6c 

10 

35 

75 

00 

0095 

D-7a 

10 

35 

85 

25 

0043 

48 

D-7b 

10 

35 

80 

25 

0084 

4.1 

D-7c 

10 

3.5 

75 

2.5 

0.047 

3.7 

D  8a 

10 

35 

85 

2.5 

0  092 

39 

D-8b 

10 

35 

80 

2.5 

0035 

D-8c 

10 

35 

75 

25 

0069 

APPENDIX  A 


LISTING  OF  THE  COMPUTER  PROGRAM  BRNMSH 


iVO  *** 
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PROGRA"  PRNMSH 

A<  I NBUT, OUTPUT » PUNCH, TAPE5S INPUT  . T APE 6s0 UTPUT .TAPE TsPUNCH > 
C  PROGRAM  ;RN"SM 

C  A< IN®UT, OUT® UT, PUNCH tT A  PE  5  s ISOUrtTA PE6S0UTP0T, TAPE 7=PONCH» 

COMMENT... 

COMMENT. ..THIS  ROJTINE  DEVELOPfS  THE  INPUT  FOR  eLE  $H3 

COMMENT...  STATUS - OPERATIONAL - IRMAT79 

CO“MENT ... 

DIMENSION  X  C  90 ) •  7(90),  AA(90>«  °9(90),  ’Hl(90)«  THFI90) 
DIMENSION  NN<  90  > ,  \P(90) 

DIMENSION  I  END  f3 ) 

OATA  LUl/5/.  LUO/6/,  LUA/7/ 

OATA  IENO ( I  )/lHE /,  IEND(2)/1hN/,  IEN0(3>/XH0/ 

DATA  PI /3. 1 4159265/ ,  PID/1P0.0/,  OEGRAO/1 .TA532925E-02/ 
OATA  2ER0/0.0/ 

DATA  IZERO/O/,  NC/90/ 

OATA  1TITLE/10HGRA1N  MESH/ 

6111  FORMAT(AIO) 

WRITE  (LUA, 61111  ITITLC 
5211  FORMAT ( 7E 1 0 .2  > 

READ  <  LU 1 , 52 1 1 )  XN,  RC,  W,  PHIDt  YA,  ALPHAD,  XL 
5311  FORMAKF10.2) 

READ  <L'JI,5311)  BURN 

PHIsPHICOEGRAD 

ALPHAsALPHAG*OEGRAO 

THETAOsPID/XN 

THETAsTHETAD»OEGRAD 

IMAXS1A 
JMAXS5 
NNC=39 
IBCs  J 
SCALES] .0 
NLAYS5 
Z=2ER0 
DZ=2.0 
ZmaxsS.O 
DO  111  Isl.NC 
X( IlsZERO 
Y<I)=ZER0 
NN(I»sIZERC 
NP< 1 >  =  I  ZERO 
A A ( I ) rZER  C 
BBIIIsZERO 
THICIIsZERO 
THFtllsZERO 
111  CONTINUE 
EPSLONsO.l 
RCWYAsRC-W-YA 
PINPHIsTHETA-PHI 
SFPHlsSINCPINPHIl 
CFPhIsCOS(PINPHI> 

CELT A*P I *0.5* THETA-PHI -ALPHA 
SFDELsSIN(OELTA) 

CPDELsCOSTOELTAI 

SPHlsStNIPHl) 

CPMIsCOSTPHI) 

SFAPsS INI  PHI* ALPHA) 

CFAPsC OS (PHI* ALPHA) 

BET A«PI« 0.5-PHI- ALPHA 
SBET  AsS IN( PET  A > 

CBET  AsCOS ( SETA ) 

APTsALPHA«0.5*PHI»0.5*THETA 
SAPTsS IN  I  APT) 

C APTsCOS  <  APT) 


15 


not 


YC=<RCWYA»SFPHI-YA*CF0EL-XL»SF2EL»/SIN<  ALPHA *PH I- TM£ T A I 
YC8=YC-SURN 

rSTAR*RCHYA*SFOHl/SIN<  ALPHA »PH I -THlTA )-YA 

Y9  =  Y  A*SUR  N 

B2=BURN-YC 

93=BURN-YSTAR 

hs(YSTAR-YC-92)*SIN< ALPHA-TmETA>/SFDEL 
m1=<YSTAR-yC-B2>oSIN<4LPHA*PHI-ThETAI/SFDEL 
I  FTPmI.lt. theta. ANC.BURN.LT.YC)  I8C  =  7 
*211  FORMAT <  *  1 5«  F10.2.I5) 

WRITE  (LUA.6211)  IHAX,  JMAX,  NNC«  I“C ,  SCALE*  NLA Y 

201  CONTINUE 
1  =  1 

NN  < I >  =  1 
NPCI»=THETAD 
A A ( I I =EPSLON 
09(1 ) =E°SLON 
THF ( I ) -THET  AC 

202  CONTINUE 
1  =  3 

X< I»=RC 
NN( 11=2 

203  CONTINUE 
I  =  A 

NN  <11  =  3 
N°(I)=THETAC 
A A ( I ) SnC 
BB(I)=RC 
THC ( I >  =  THET  AO 

204  CONTINUE 

1=6 

X< I>=RC 

Y(I)rR-«TAN(THETA) 

NN( I >=4 

205  CONTINUE 

1  =  7 

NNU  )  =  0 
NP(I)=THETAO 
AA(I»=RC-M 
9«l]»tll]l 
THF ( I l=THETAO 

206  CONTINUE 

1=6 

CC=(RCMYA*SPHI*YP»SFAP*XL»S9ETA)-(RCWYA •sphi*yb»sfap> 
00=IRCWYA«CPHI*YB»CFAP)-(RC JYA »CPH I ♦ YB» CFAP-XL*CBET A> 
SLOPE l=CC/DO 

X<I)=RCWYA.CPHI*YB»CFAP-XL*CRETA-0.1 
Y(  i>=rcwya«sphi*yb*sfap*xl«speta*o.i»slopei 

1  =  9 

X(I»=RCWYA«CPHl*YB«CFAP*0.1 
Y< I)=RCWYA«SPmI*Y0»SFAP-O.1«SLOPE1 
NN(D=e 
2 0 T  CONTINUE 
1=10 

X( I )=RCWYA*CPHI 
Y  < I l=RCW*A»SPMI 
NN<  t)  =  T 

np ( I >* alpha; 

AAU)=*‘ 
asmsYB 
ThI ( I IsPHIO 
THF(IlsPMlO»ALPHAD 
20R  CONTINUE 
1  =  12 

G5S=RCWYA«SPHI*YB«SFAP»XL*SAETA 

HHHsRCWYA»CPHI*Y9*CFAP-XL»CRETA 
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SL0PC2*0GG/NMM 

xit  >*ncwy a«cpmi»yp«cfap-*i«c9ETa*o.i 
YU  >**CUYA« SPH|*YP«SFAP*Xl» 4»ET 4*0.1 «S19PE2 

nnump 

20q  CONTI NUE 
1  =  13 

*U)=4C4YA.CPMI*YP«CFAP-»L*CEETA*YCB«SeETA 
T(  I»*4C4YA.SPMI*Y9«SFAP*XI..?9CTA*YCB«CP:TA 
NNU>*9 

NPU I=AIPHA3*PhI0-TmCTAD 

AAUlaVCP 

B3UHYCP 

THJ(  I )*«I0*TMETA0-1.0 
THFU»*PID*PHI0*4I.omA0M.0 
21 C  CONTINUE 
1*15 
mime 

YU)*BC«TAN<PHI) 

NN  U  >  =  1 0 

211  CONTINUE 
1*16 

NN( 11*11 
NPU  »=TnETA9 

M(I  »*URCWYA*CEPMI*Y?*<!fDEL  » • • 2 ♦ < *CM YA *SFPMI -Y 8 ECFDCl > • • 2 > * • 0 

BBU  1:111 1 1 

TMFUIsTMCTAD 

212  CONTINUE 
1=1? 

KIlltRC 

YU>*RC«TAN(TMETA«0.5» 

NN ( 1 1  *  1  2 

215  CONTINUE 
1*19 

NN< 1 1  *  1 3 
NP(I )=TMETA0 

CCC*(PCUY»«CFPHl*TP»SF0EL-YL»CF0CL>«»2 

000=MCWY*«SFPM1-YB»CFOEI.-XI«SFOEL1««2 

AAU»*ICCC*000>*«>0.« 

BBIIXUII) 

TMFU  IsTHETXO 
214  CONTINUE 
I  *2C 

EE  =  XAU9|»SINUHET»*0.5>>*M16)*SPHl 
FF*AAU6l«CPHl-AA<l<»>*COSUHCTA*0.5> 

SL0PE3*CE/FF 

XU>  =  AAU9>*C0S(THET»O.5>-0.1 
YU  »«AA(19>*S1 N< THETA* 0.5 > ♦O. 1 • SLOPE 3 
1*21 

XU):AAI16>*CPH1*0.1 
YU»*AAU6)«SPPl-0.1«St.0PF  3 
NNU  >  =  14 
21*  CONTINUE 
I  *22 

NN<  I  >*19 
NPU»*ThEYAC 

AA(  n*EPSLON»IAA(l<l)*AAU>)*O.T* 

BBUXAAU) 

THFIDsTHETAO 

216  CONTINUE 
1*23 

GG*AAU6)*SPHI-AA|T)*SIN(PHI*0.5> 

HH*AAU)*COSIPHI*0.5>*AA(16)*CPhI 

S10PC4*GG/MM 

X<  I>*AAU6>*CPMt-0.  t 

YU)*»AU4>*SPMl*0.1«9l.OPE4 

1*24 
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xiiizaaiti. cost  phi* o.si.q.i 

Tt I)*AA( 7»«SIN(PhI«0.5I-0. 1»S10PC« 

NN(I)*16 
217  CONTINUE 
1*25 

NN(I>=17 
NPtI )*TMCTAD 

AA(I)=AA(19)*(AA(16)-A 1(19)1*0. 75 

eeiiiiitu) 

THP  C I  I  =  THE  T  AO 
21 P  CONTINUE 
1=27 
X(  1  )  =  *C 

rt 1 )sRC*T*N(PHl«0.5l 
NNtl)*l* 

219  CONTINUE 
22C  CONTINUE 
1  =  28 

VV*RC»(YB.SIN< ALPHA  I /(SCyYA.YO. COS (ALPHA) 1 1 -RC«Y»«SPHI ♦YB«SF »P 

lMI*NC-ltC2YA.CPHl*Y?»CCAP 

SLOPE5=VV'WU 

tt I)*9C«Y»*CPHI»TB«CF*3-0.1 
Yt  I)«ltCWYA«SBHl»vp.$FAP-0.1*SLOPC5 
1*2* 

XtllMC 

Yt I)*NC.tY«.SINt ALPHA) /t9CUYA*Y».C0S( ALPHA) 1) 

NNt 11=20 

221  CONTINUE 

222  CONTINUC 
1*51 

X ( 1 1 «NC 

Yt I>=(AC.tY3.SINtAL“HAI/tAC«YA.Ye*C0$(ALPHA> )) )»0.5 
NNt 11*22 
225  CQNT lNUr 
22*  CONTINUC 
1*52 

UUU  =  Y (51 )**2/(8C**“) 

NNW*(K(10>*t( V*».2).(1.0«UUU>-( M 1 0 1 • «2 1 »UUU) • • 0 .5 1 / ( 1 . 0 *UUU 1 

litiMiOttli/IC 

YTT*»«tl9).StNttNeT».«»91*?2i 

*W**WNtt-AAt  l*)»COS (THETA* 0.f> 

SL0PE6=yYY/XXX 

X(I)=AA(19>*COS(THCTA.a.S)-C.l 

Yt I)=AA(19).SIN( THETA. 0.5) ♦0.1*SLOPC6 

1*55 

X(I)=UUU»0.1 
Y(I)=2Z2-0.1«SLOPfS 
NNt 11=2* 

225  CONTINUC 

226  CONTINUE 
1*55 

SlOPET=(Y(15)-YC««SAPT>/tXtl3)-YC°»CAPT  1 
X(I)=x(15)-YCB*CAPT.0.1 
Yt I)=Y(13)-YCB.SAPT*0.1«SLOPE7 
NNt  11=26 

227  CONTINUC 
22  *  CONTINUC 

1  =  56 

XtI)=Xtl5)-YCB*CAPT-0.1 
YII)=Y(15)-TC?»S APT-0. 1«TAN (THE TA.O. 25) 

1*57 

X( 1 1 =AA( 16) 

YtI)=Y(13)-YCd.SAPT.(AA(16)-Xtl5>»YC«.C*PT).TANtTHETA.0.25) 

NNt 1 1  =  28 
229  CONTINUE 
250  CONTINUE 
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1  =  50 

AD=XI15)-»C-*CAPT 

AC=YU3)-YC!<*S43T 

kfz  i»  (  1C.  ) 

A-,=  T»NITmCTA*0.25> 

AXs4G**2. 0*1.0 
*»i2.0*'AD**C**r)> 

»/s40**2.0*«t**J.O-»c**?.0 

By=AW*4(, 

CU=4AI  T)*Sl*!(T(*ETA*0.5>-(  Yl  l  31-YC®*SAPT  *BW) 
DU  =  A*IT>*COS(TmETA.''.')-(X(  1 3)-YC?*CAPT  *AM> 
siOPER=Cu'&y 

X|l>  =  XU3)'YC**CAPT*AU-0.l 

Y<!)=YC13)-YCa*SA®T.3w-j.l.SL0Pr4 

I-co 

X<I»s4A<7j.COS(THfT4.a.b»*0.1 
Y I 1 IxAAt 7>*S1NC’HET**0.5>*0.1 *SLO®E« 

NN  < 17-30 
ill  CONTINUE 

?s:  continue 

I=4n 

xii>=rc*o.2k 

Yl  I  >sRC*0.25*TANU«!:T4.o.M 

1^41 

X<l)=RC 

Y  I  I  >rRC*TAMTMET  4*3.5) 

NNl I )=32 
235  CONTINUE 
2  34  CONTINUE 
1x43 
xi  disc 

Y I  1 )  sR  C*T  A\ I T  HET  A  «  0 . 35 ) 

NNl 1)554 
235  COMlNOr 
2bb  CONTINUE 
1  =  45 
X( t )  =  RC 

Y(I)sRC*T44|(ThET4*0.58) 

NN I  I >  =  ' b 
2  57  CONTINUE 
25s  CONTINUE 
1547 

X  1  1  )  =  4  C 

Y I 1 >  r 4  C*  T  A\  I  T  MET  4*3. 75) 

NNl i)*?a 
2  3  4  CONTINUE 

240  CONTINUC 
I  34  R 

AB=AAIlb)*SINITMCTA*0.5)-IY<l J)-YC?*SAPT) 

AC  =  4A114>*CnSCTHETA*0.4)-««U3>-YCb*CAPT) 
SLCPERsAB/AC 

XI  I)=XI13>-YC°*CAPT-0.1 
YII)=Y<13>“YC5*SAPT-0.1« SLOPE® 

1344 

XU)sAA(lb>*COS<THCTA*0.S>*0.1 

YU)  =  AAU5)*SIN«THET4.  0.5)*0.  l*SLOPrN 

NNI1)=40 

241  CONTINUC 

242  CONTINUE 
1=50 

XII)*RCWYA*IYA*YC*B2> 

1=51 

XI I) =RCWYA*<Y4» YC*3?)*r OS I ALPHA/ 3* 0) 

Tills <YA*y:*92>*SIN<AI°m»/J,0) 

1=52 


XU»*XC«x»»IXA»XCy.’l«C'l$l»t.P**A*J,0/S.C> 

Ttl»><Y**YCM;i«$lXIM.PH*«t.O/J.C> 

I  *?  5 

X|||t»C*X»»IX»*XC«!>‘,l*COSt»lP'«»» 

TII>s|x»»YC«P?l*StNIAl»HM 

!»!>• 

XII)*XI?>JI»»**C0$IPl*0«b**LX,**l/J»0 
*l|»*YIX5l»*«SIN|Pt*0.»-Al.P‘*»>/.,.0 
I  «*5 

NNI t  »*♦? 

2*J  CONTIM"- 
I*3fc 
*  (  I 

NN I  1  I ** ' 

Np I  I  )  * *LPH*P 
AAI  I  )*Y»«YC»«l  J 

mi  n  *»*<  i  > 

ThF I  I >s»LpM»0 
CON*INur 
I  «*  r 

XI 

tsr  S 

XI  I >*XIS1 >/l. 125 

YII»«YIM»M.1J? 

xi n  *x«52»/i. 12* 
yi i»*yi5»»/i.»j« 

I  *b0 

xi  IHXI5SI/I.U? 

y  1 1  i>tij3im<u; 

I--M 

«in>«IJMM.l2! 

TiDniMm.  u« 

1«6? 

x<l>*xisS»/WI25 
yuxyiM)Mam 
NNI  t  >»*• 

2*5  COMUHUC 
1*6? 

xi  n*scxYA 

NN  C I»*X5 

*P< l*=Xl°M»3 

MMIXtlSO/MJ* 

MIIHUIII 
TNM  ICILPHtC1 
2X6  CONTINU*' 

I  -fc« 

XIII  *XCyYA»IYA»YCx«2* 

1*6? 

xi  n  *uc-y*YCxp*»«cp*<i 
y«  i»*i*c-w*yc«P2i  •$«**! 

l  *‘6 

xil»*xCwY»*CP-«l»IY»«YC«,2>«:C'xi*>Mt»»l*MA23.0> 

y(l»*xcyy»«s»«i*«y»*TC»P?>«s!HtP**i*iv.pyu<3.o» 

!*‘T 

xi  i > *xr yTt«cPHl* (ywy:»l  ;i« r'Si0*-!* h.ph 

yU)*xcxy*«SPHj*iyi»yf»?>»«xtNipyt.*tPH*«f.o/>.o» 

l*AX 

xi  n*xc»y*«,‘PMj*iy*»y»P?»«cc*f 
yil)«xcxy*»s»Hi»iy*»yr«x?>«xr»p 
t 

XI  I)*X|kM-Ht«CXtT*/I.O 

Yl  l  >*YIMI*m«SXtTAFf  .0 

t  *  ?  o 
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X(Il=X(6R>*Ml,CeETA 
VI  I  »  =  Y  (  6>fl  )  «H  1  .C9ET  1 

MN<  I>**5 
2*7  CONTINUE 
1  =  71 

NNU  »  =  *  7 

NP(  I  >=PhIC 

AA  I  I  »sSC-W*YC«r  3 

se I  I t : A  A I  I  I 
TMFIDsPHlC 
2*8  CONTINUE 
TsT  r 

«<  I)  =  X(M>/1.125 
1=73 

Xlt)s«(65>/1.125 
Y«I»=Y(b5>/1.125 
1  =  7* 

K(II=X<b6l/1.125 
Y< I)=Y(S5)/1.125 
1=7* 

Y I I)=X(67>/1.125 
Y( I >=Y(67>/1. 125 
1  =  75 

X(I)=X(681/1.125 
Y(I>=Y(68»/1.125 
1  =  77 

»II»=»I&9>/1.125 
Y(I)=YI69>/1.125 
I  =  7  8 

XI  l)=XC70>/1.125 
Y( 1)=Y«70)/1.125 
NNCI »=*6 
2*«  CONTINUE 
1  =  79 

XI Il=RCWYA*CPHl 
YII)=RCWYA»SPHI 
NNI I >  =  A  9 
NP(I >=ALPH«D 
AAI I >=YA*YC»B3 
33 ( I >  =  AA(  I  ) 

THF  1 1 >  =  ALPHAO 

250  CONTINUE 

251  CONTINUE 
I  =80 

NNI  I>=51 
NP 1 1 ) =PH 1 0 
AAI I)=AA(71 1/1.125 
63(1 )=A»< I ) 

THFC IlsPPIP 

252  CONTINUE 

253  CONTINUE 
I  =  31 

XI I»=RCWYA*CPHI 

YII)=RCUYA«SPHI 

NNI I)=53 

NP I  1 1 = ALPHA 0 

AA(1)=AA(7«)/1.125 

BPII>=AA(I) 

THF 1 1 )  =  ALPH AD 
299  CONTINUE 

00  331  1=1, »1 

WRITE  (LUA.6331)  X ( I > , Y 1 1 ) , NN I  I ) ,NP I  I  I , AA ( I > , PR ( I ) , THI ( I > , THF I I» 
331  CONTINUE 

F  331  FORHAT(2F10.*,2I5,*F10.A* 

INl=-l 
I I=2ER0 
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I  J  =  l 
I«  =  2 

II  =  ! 

IH  =  A 
1N  =  R 
IP  =  G 

10  =  7 
IM« 

It;: 

01  =  13 
JJ=  1 1 
JR  =  1 2 
Jl-15 
JN:l» 

JN  =  15 
Js  =  lf 
J0  =  17 
JR  =  1? 

JS=1  = 

•<1=23 
KJrJl 
*<  =  ?2 
■a  =  23 

«P  =  2A 
*N  =  25 
*P=t6 
R3  =  ?7 
*R  =  29 
<«sJ9 

1 1  =  3  C 
t  J=3  1 
L<  =  32 
LL  =  '3 
lp=’a 
I.N='3) 

LP=’G 
1.0  =  37 
LR  =  3F 
LSs'R 
l»I  =  AC 
*J  =  A1 
PK  =  A? 

NL  =  A  J 
WRr  AA 

<*N  =  «5 
>'P  =  «6 
•  Q:t  T 
«R  =  AR 
KSsaR 
N 1  =5  0 
N  J=5  1 
NR  =  52 
NL  =  R! 

I  73  jc  =  r 

lP(eUR!'i.Lt.rC)  1 20N?  =  1 

IF<=UR\.GE.»C  .ANO.  SuRN *L T  .  Y S T AR  )  I20NE*2 
IM-UR»  .GE. TATAR  .ANO.  RUR'.LE.W)  I’ONfrS 

ici  uone.nc.o  50  to 

A'-P  C  3\  T  I  NUf 

WRITE  ( LUO«  6A  3  p>  sURM,rC*Y°TAR ,y 
b*3S  FORMAT  mi  20 NT.  ERROR  .  /  t  A  U  X  .  '  T  0  .  A  >  > 

&o  TO  PP9 
•*0  CONTINUE 

IF  (PHI.LT.THCTA)  50  TC  A  A  ? 

6373  FORPAT  I5IR) 

WRITE  U.UA,  63731  I  J  .  I  J .  1  J  «  1  '  .  I  N 
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**  1 
b*.  ' 


s  jj; 


*♦3 
b*.  1 


RRE 

SR.2 


WRITE  luU.bJ’J) 
a»l'f 

CONT  !Njr 


F  OR  “A 

’  (1*1‘> 

WR  I  Tf 

(LUA.bR.S) 

W  R  I  T  l 

<  LU*  ,b»K * ) 

WRITE 

IE  U  A  «  b  *  R  3  ) 

■  RITE 

(LuA.h*. 0) 

WR  ITE 

(  ,.i  A  .  h  *  b  5  ) 

WRITE 

ILUA.bR. 3) 

WRITE 

(LUA.bR.3l 

■  RITE 

(LUA.bR.3) 

WR  I  ’  E 

«  LU*«  SR. 3) 

WRITE 

(LUA.bR.3) 

WRITE 

(LuA.h*. 3) 

WRITE 

>L JA.bR.3) 

WRITE 

(LU»«b.*3) 

SO  TO 

*? : 

COM  I  NU  r 


IF ( I20NE.NE  .1  )  '■.0 
*ORRAT ( ; 15) 

WRITE  (LUA,b332> 
WRITE  (LU).b33J) 
WRITE  <LUA,bSS2) 
WRITE  <LUA,b3S2> 
WRITE  ( LUA  «  N332 ) 
WRITE  (LU*.b332> 
WRITE  (LUA,b332) 

C  ONT I  NU  r 

Ic(PmI,EO.O.O)  >0 
FORMAT ( 1  A  15) 

WRITE  (LUA.6RR1) 
WRITE  (LUA.bRRl) 
WRITE  (LUA.6R.1) 
WRITE  t LUA.bRRl) 
WRITE  (LUA.bRRl) 
WRITE  (LUA.Frr 1 ) 
WRITE  ( LU  A , bRR  1 ) 
WRITE  (LUA.bRRl) 
WRITE  (LUA.bRRl > 
WRITE  (LUA.bRRl) 
WRITE  (LUA.bRRl) 
WRITE  (LUA.6RR1) 
WRITE  (LUA.bRRl) 
WRITE  (LUA.bRRl) 
WRITE  (LUA.bRRl) 
WRITE  (LUA.bRRl) 
WRITE  (LUA.6RR1) 
WRITE  (LUA.bRRl) 
WRITE  (LUA.b**l> 
00  TO  R51' 

CONT INUr 
FORMAT  (IRI5) 
WRITE  (LUA.bRRi) 
WRITE  (LUA.bR.D 
WRITE  (LUA.bR.D 
WRITE  ( LU  A , b. R  ? ) 
WRITE  (LUA»bRR2) 
WRITE  (LUA.bR.D 
WRITE  (LUA,b**2) 
WRITE  (LUA,bRR2> 
WRITE  (LUA,b*R2) 
WRITE  (LUA.6RR2) 
WRITE  (LUA.bR.2) 
WRITE  (LU*. b*R2> 
WRITE  (lUA.bRR2) 


IK.JI.IJ.JI.IN 
IL  .  JJ.  I  J.  J'’.  I  N 

IJ.lJ.Iv1  .IS.  Ik 


IR.JJ.LR.IJ.lJ.l 


iK.l J.lS.lS.lL.LN.JJ.LP.l J.lJ.l 
Il.U.IL.lS.I-.LP.JJ.LR.IJ.lJ.I 
I  -  .  I  J  .  I  N,  I  E ,  I  \  .  L  R  •  J  J  .  I  I  J  .  I  J  .  1 
1\'.1S.1J.J1.1k,1K.1N,LN,JJ,1K,1 
!•>.!?,  IK.J1.IL.  I", IN.  IP.  JJ.IK.I 
I".  .IS,IL.J1.IR.LP.IN.LR.JJ.IR.1 
Ih. IS, ]N.JI. IN. LR.IN.1“. JJ.IK.I 
1 S. JI , I J.J«, IK, IK.IL .L*. IN.IL.I 

JI, J!.IK,J“,IL.L*.IL.LP.IN.IL,I 

JJ. JI.IL.JN.I", LP.IL.LR.IR.IL.  I 
UK.JI.IN.jm.IN.LR.IL.IN.IN.IL.I 
IRl.II.IItll.II.lI.lI.lI.II.il. 


TO  »>« 

IJ.IJ.IJ.1n, I* 

I J. IN.l  J.I'.lK 

IK. IA.IL.IN.IN 
IK.IP.IL.IR.1L 

IK. JI.IJ.JI.IL 

IL .  IP. 1 N, J«.  IN 
IL . JJ. I J. J". IL 

TO  R*5 

IJ. IJ.IJ.IR .IK.IK.JL.JK.IJ.IJ.I 

IK. IJ.IK.IN.IL.JK.JL.IR.IJ.IU.I 

IL. IJ,IL,IM«IM.IR*JN.KP»IJ.IJ.I 
IK«IJ.IM.l“.lR»KP,JN,lv«lJ,lJ»l 
IN.IN.IJ.1S»IK.1K.JJ,JW,JL»IJ.1 
IP.IR.IL.IK'.IR.IR.IS.KP.JN.IK.I 
l 1. IN, IN, IN* IN.KP, IS.I R.JS.lK.I 
IR.IN.IK.IF.IL.JN.JJ. IP. JL.IK.I 

IR. IS.IJ.JI.IK.IK.1N.jp, JJ.IK.I 

JI.  IS. IK, JI.1L.JP. IN, 1G, JJ.IK.I 

JJ . IN. IL.IS.IR.IP.JJ.KR.IS.IL. I 

JK.  IN, IN. I R, IN.KR. JJ.I «, IS.1L. I 

JL . IS.IL. JI .1*. 18. IN, L 1. JJ.IL t I 
JR. IR. IN. JI.IN.LI.IN.1R.JJ.il. I 

JN,  JI. I J, JR, IK.IK.IL, JR.IN.il,  I 
JP. JI . IK, JN, u . JR . IL .J  1, IN.IL ,  1 

JO. JI.IL, JR, M, Jl.lL.LK.IN.IL, I 
JR. JI. I*. JR. IN.LK.IL, IN, IN.IL,  I 
tRl.l  1.1 1  ,11  .1 1.1 1  .11.  II  .1  1  »I  I. 


lJ.IJ.IJ.IN.lK.lK.JL.JK.IJ.IJ.I 

IK. IJ.IK.IN.IL.JK.JL.IR.IJ.IJ.I 

IL. IJ.IL,I“,I*.1R»JN»RP,IJ,IJ,I 

IR.IJ.1“,I",IN.KP,JN,1R,IJ,IJ,I 

IN. IN.IJ.IS.IK.1K, JO, KR.JL.IJ, I 
IP.IR.IL.IV. IR.IR.IS.KF.JN.IK.I 

IO. IR,IN,IN,1V«KP,IS,I“,JN,IK,I 
IR. IN, IK.IR.IL, KM, JQ, IP, JL, IK, 1 

IR. IS.IJ.JI.IK.IK.IQ.KR.JO.IK.I 

JI. IS.IK.JI.IL.KN.IQ.IF.JQ.IK.I 

JJ.  IN.IL. is. I v. IP, JQ.NI.IS.IL,  l 

JK. IN.IR.IR.IN.RI »JQ»I ».IS. IL.l 
jl.is.il»ji.i«»ip.jj»ri»js.il»i 


I, 


I. 
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W*»  I  T  t  (III*.  *>**.')  J«.  I  $  .  I**,  Jl  .  IN.M  I  •  JJ.  1  M.  JO.  IL  .  I  1  •  I  I  •  I  I  t  I  l 
WRIH  T  LU  A , bA  A  0 )  JN,Jl.IJ.JM.lK,IK,lL,KK,|3,lL,ll,tl,IUII 
HHIH  (  LUA.bA*.')  JR,  JI,  IK.  J“,  II, KK,  IL.K  1,  10. IL,  I  I,  I  1,1  1,11 
J»IU  (LUA.bAAI’T  JO.  Jl  ,  1 1  ,  J*,  1  k.k  |  ,  IL  ,L  •.  JJ.  I  L,  1 1  ,  I  I  ,  II  .  1  I 
WHITE  TLUA.fAA})  Jl, Jl ,|M, Jl, IN. IK, | L . | M, JJ.IL, II ,1 1,1  I ,1  I 
WRIU  llU*,M*r>  Ml.ll.ll, I], II. ll.II.Il.il, 11,11, II, II, II 
GO  TO  *10 
»AA  CONUNU! 

I»  (  IZONC.NT  .?>  GO  TO  AAb 
633A  FORMAT  (-'ll)) 

WK|TF  (LUA.633A)  IJ.IJ.IJ.JJ.I* 

WRITE  UUA.fS.'A)  IK  ,  JK,  |  J,  JK,  IV 
WRITE  (LUA,b33A)  IL.JL.IJ.JM.IN 
Ml  CONTINUE 

iMPHi.ro. o.o>  r.o  to  a  a  11 

MAT  FORMAT  (1*1  til 

WRITE  (LIlA.bAATT  I  J  ,  I  J  ,  I  J  ,  J  J  .  I  N  .  I  A  ,MH  ,  I  M,  I  J  .  I  J  ,  1 1  ,  I  I  ,  1  I  ,  1  J 
WRITE  (LUA.hAAT)  I  a  ,  JJ,  I  J.  JK,  IN,  IK.mp,  I  «,KA  ,  |K,  I  I  •  I  I  ,  I  I  ,  I  I  . 
WRITE  (LUA.MAT)  IL.JK.1J,  JM,|N, IN, R.U'. "P»IL, 11,1  1,11. 1J 
WRITE  (LUA.hAAT)  I  R|  .  |  |  ,  l  |  ,  I  |  ,  I  I  ,  I  I  ,  |  I,  I  I,  I  1,  I  I,  |  1,  I  1,  I  I  ,  I  | 
GO  TO  A  0  0 
AA'l  CONTINUE 
hAAN  FORMAT  (|A|N) 

WRITE  (LOA.MAM  I  J  ,  I  J  •  I  J  •  J  J,  I  N  ,  I  •  ,  MM  ,  I  M,  I  J  ,  |  J  ,  I  1  ,  |  |  ,  1  I  ,  I  J 
WRITE  (LUA.hAAM  I  K  .  JJ,  I  J,  JK,  |  N,  IK.MK,  I  M.MM,  IK,  I  I  ,  I  I  .  I  I.  I  I 
WHITT  (LUA.hAAN)  1 1  ,JH , I J, JM, IN. I K, IL ,1 M,MK, IL, 1  I , 1  I . I  I , I J 
WRITE  (LUA.hAAN)  I  *  I  ,  I  1  .  II  ,  I  I  .  I  I  ,  I  I  ,  I  1  ,  I  I  ,  I  I  ,  I  I  ,  1  1  ,  I  I  .  I  1  ,  II 

GO  TO  A  Ml 

AAb  CONT INUr 

ir  (  l ZONE  •  N!  .  S *  NO  TO  A  IP 
b  3  3  b  FORMAT  (NIK) 

WRITE  (LUA,b33b)  1 J . I J.l J, JK, IN 
WRI’C  ( L  U  A , b  3  3b  >  IK , JL . I J. JL, IN 
WRITE  (LUA,bS.3h>  IL.JM,  I  J,J**,]M 
AM  CONT  I  NUr 

lFEPHl.LU.0.0)  NO  TO  A N 3 
b  A  A  H  FORMAT  (IAIN) 

WRITE  (LUA.bAAH)  | J, | J, I J, JK , IK, I K ,NJ,J  | , I J,  I j, I  1 • I  1 , I  I , I J 
WRITE  (LUA.bAAH)  |K,|J,|K,JK,1N,JUNL,1M,IJ«|J,II,II,1UIJ 
WRITE  (LUA.bAAH)  I L , JK , I J, JL , 1 K , I K , N 3 , J  I  , N J , I K , I  I, I  I  • I  I, I  I 
WRITC  (LUA.bAAE)  I M  ,  JK  .  I  K  ,  JL  .  I N.  J  I  ,MR.  I  M.NL  ,  l  K  ,  1 1  ,  I  1  .  1  I  »  I  I 
WRITE  (LUA.bAAH)  I  N.  JL  ,  I  J,  JM,  IK,  IK  •  IL  ,  J  I  .Ml),  IL,  1 1  ,  I  1  ,  I  I  ,  I  1 
WRITE  (LUA.bAAH)  |P, JL. IK, JK.IN.JI.IL.IM. MS, IL.II.II.II.il 
WRITE  (LUA.bAAH)  I M | , 1  I , 1  I , 11 , I  1 , I  1  .  I  I  ,  1 1 , I  1  ,  I  1  ,  II  .  1  I  .  I  I  ,  I  I 
30  TO  AN(1 
AK3  CONTINUl 
bAAb  FORMAT  11,13) 

WRITE  (LUA.bAAb)  I J, | J, 1 J, JK. IN, I K , MN , I “, I J, I J. I  1 . 11 , IT , 1 J 
WRITE  ( LUA.bAAb)  I N , JK , I J , JL , I N , I K • ML , 1 M.MN , 1 K • I  I , 1 1 , I  I ,  I  1 
WAIIf  (LUA.bAAb)  IL • JL, I J, J A, IN. IK.IL.1 m.ml.IL,! UI I, 1  I, I  I 
WRITE  (LUA.bAAb)  Ml ,i |, | |, | |, | | , | 1 , | |, U ,11 ,  I  1.1  I  ,  1 1  ,  1 1  ,  1  I 
ANC  CONTINUE 
b  A  A  A  FOR  TAT  (MO. A) 

WR 1 TC  (LUA ,h*A A  >  ’ 

7-ZaOZ 

IT  (■•.&!. .’MAX)  GO  TO  >«'> 

ANN  CONTINII' 

1 F 1 1  ZONE .EG. 1 )  30  T  T  *  A  3 
ANF  CONUNU'’ 

IF (  IZONC.IO.?)  -0  TO  AA  » 

•  NT  CONTINUl. 

1 F  <  I  ZONE . E  0 . 3 )  GO  TO  «N| 

ANA  CONTINUE 

IFEPHI.LT. TmETA)  GO  TO  a  A l 
b 33 3  FORMAT ( 3 A  l  ) 

WRITE  (LUA«b333)  ( I E NO ( I ) , 1 3  I , ' ) 
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APPENDIX  B 


EXAMPLE  THREE-DIMENSIONAL  FINITE  ELEMENT  MESH 
GENERATION  INPUT  FOR  STAR  DESIGN 


I8ECKDIMG  PAGE  BUMK 


GRAIN  NESH 


Zone  I  Burning 


14  5 

39  7 

1.00 

5 

0.0000 

0.0000 

1 

25 

.1000 

.1000 

0.0000 

25.7143 

_ n-nnnn 

.-  0.0000 

0 

0 

o.oaoo 

0.0000 

0.0000  . 

8*0  000. 

10.0000 

0.0003 

2 

0 

0.0000 

0  .0000 

0.0000 

0.0000 

- 0*0000 

0.0000 

3 

25 

10.0000 

10  .0000 

0.0000 

25.7143 

0.0000 

0.0000 

C 

0 

0.0000 

0.0000 

r.oooo 

0.0000 

10.0000 

4.8157 

4 

0 

0.0000 

0.0003 

0.0000 

o.oooc 

0.0000 

0.3000 

c 

25 

6.5000 

6  .5000 

0.0000 

25.7143 

- 2*8348 

-  1*4 122 

0 

0 

0.0000 

0.0000 

0.0000 

0*0008 _ _ 

6.2424 

.5  045 

6 

0 

0.0000 

0  .0000 

0.0000 

0.0000 

_  6.0000 

0.0000 

7 

75 

.5500 

.5500 

0.0000 

75.0008 

0.0000 

0.0000 

0 

0 

0.0000 

0  .0000 

0.0009 

0.0000 

-  _  3.054H 

1.4322 

a 

0 

0.0000 

0  .0000 

0.0000 

O.OCQO- 

2.9664 

1.4285 

9 

49 

.0447 

.0447 

204.7143 

256.0000 

- 0.0000 

-0*0000 

0 

a 

o.oooo 

0 .0000 

0*0080 _ 

n.ocno  .. 

10.0000 

0.0000 

10 

0 

0.3000 

0  .coco 

0.0000 

0.0000 

_ 0*0000 

0.0000 

11 

25 

6.1653 

6.1653 

O.OOOO,. 

-25.714! 

C.0000 

0.0000 

0 

0 

0.0000 

0  .0000 

0.0000 

0.0000 

_ .10*000.0 

2.2824 

12 

0 

0.0000 

0  .0000 

0.0000 

O.OQQO 

0.0000 

0.0000 

13 

25 

3.2634 

3.2634 

0.0000 

25.7143 

- 1.0016 

—  .*7503 

C 

a 

fl.CQOU 

o  .(umo 

o*oaoo  _ 

.  n.cnno _ 

6.2653 

-.0243 

14 

9 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

15 

25 

2.4726 

2.4726 

O.OQOO  . 

25.714A  _ 

6.0653 

0.0000 

0 

0 

0.0000 

0  .0000 

0.0000 

0.0000 

6.6nnn 

0.0000 

16 

0 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

17 

25 

5.4398 

5.4398 

0.0000 

25.7143 

_ C.nnnn 

Q^fl  QQ  Q 

.  0 

0 

0.0000 

0  .0080 

0*8  0.00. 

.  .0*0888 _ 

10.0000 

0.0000 

18 

0 

0.0000 

0 .0000 

0.0000 

0.0000 

6.0.424 

.4976 

0 

0 

0.0000 

0  .0000 

0.0000 

0.0400  .  „ 

10.0000 

.8649 

20 

0 

Q.OOOC 

0  .0000 

0.0000 

0.0000 

. JL.OO  Q  0 

0.0000 

0 

0 

0.0000 

0  .0000 

0.0000 

0.0000  _ 

10.0000 

.4325 

22 

0 

0.0000 

0  .0000 

0.0000 

o.oooo 

_ 4.n«ifc 

_ .1397  . 

0 

0.. 

0.0000 

2.0008 

.  .  0*8000 

8*0000. _ 

6.5734 

.2664 

24 

0 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

O.OQOO 

0 

0 

0.0000 

0.0000 

0.0000 

0.0000 

3.0462 

1.4358 

26 

0 

0.0000 

0.0000 

0.0000 

o.oooo 

2.8462 

1.3774 

0 

0 

0.0000 

0  .0000 

0.0000 

0.0000 

6.1653 

1.7513 

26 

c 

0.0000 

0.0000 

0.0000 

0.0000 

_ 5*8194  _ 

.  .  1.7900 

0 

0 

0.0000 

O.OOLOO 

0.0000 

o.oooc _ 

6.4370 

1.3800 

30 

0 

0.0000 

0.0000 

0.0000 

0.0000 

2.5000 

.5706 

0 

0 

0.0000 

0  .0000 

0.0000 

0.0000 

10.0000 

2.2824 

32 

0 

0.0000 

0.0000 

0.0000 

0.0000 

0,0000 

0.0000 

0 

0 

o.oooo 

0 .0000 

0.0000 

0.0000 

1C. 0000 

1.1267 

34 

0 

0.0000 

0  .0000 

0.0000 

0.0000 

-  a.onnn 

-  0*0000 

Q 

0 

0.0000 

0.0000 

0.0000  . 

.  8*0008--. 

10.0000 

2.2824 

36 

0 

0.0000 

C .0000 

0.0000 

0.0000 

Q . 0000 

0.0000 

0 

0 

0.0000 

0  .0000 

0.0000 

0.0000 

10.0000 

3.4992 

IB 

0 

0.0000 

0  .0000 

0.0000 

0.0000 

_  .  8*8462 

1.1892 

0 

0 

0.0000 

0  .0000 

0.0000 

0.0000 

6.1107 

1.3714 

40 

0 

0.0000 

0 .0000 

0.0000 

0.0000 

_  6.5500 

..  8. 0  0  40. 

0 

0 

0.0000 

0 .0000 

o.oopo. 

0.0001 

6.4985 

.2324 

0 

0 

0.0000 

C .0000 

o.oooc 

o.oooc 

6,5535 

.4213 

0 

0 

0.0000 

0 .0000 

0.0000 

0.0000 

6.1424 

.5313 

0 

0 

0.0000 

0 .0000 

0.0000 

0.0000 

.  .  *,5?35 

.9650 

0 

0 

0.0000 

0  .0000 

0.0000 

0.0000 

2.9046 

1.3988 

42 

0 

0.0000 

0.0000 

0.9000 

0.0000 

_ 6.0000 

_.  5.0000  . 

43 

.75 

-2.2899 

-2.2899 

0,0000 

75.0000  - 

5.8222 

0.0000 

0 

0 

0.0000 

0  .0003 

0.0000 

0.0000 

_  5.7764 

.2066 

0 

0 

0.0000 

0.0000 

0.0000 

0.0000 

5.6476 

.3745 

0 

0 

0.0009 

0.0000 

0.0000 

0.0000 

_ 5*9.594 

.4722 

0 

0 

0.0000 

0.019  0 

.  0.0000 

0.0000 

4.0209 

.8578 

0 

0 

0.0000 

0.0000 

0.0000 

0.0000 

2.5819 

1.2434 

44 

0 

0.0000 

0  .0000 

0.0000 

o.oooo 

29 


*• 

C030 

0.0000 

95 

79. 

-2.035* 

-2 .035* 

0.0000 

75.0000 

6.5500 

0.0000 

0 

0 

0.0000 

0.0000 

0.0000 

0.0000 

6  • 

5500 

0.0000 

0 

0 

0.0000 

0  .0000 

0.0000 

0.0000 

6.9985 

.232* 

0 

0 

0.0000 

0.0000 

o.oooo 

0.0000 

6  • 

3535 

.9213 

0 

0 

9.000C 

0.0000 

0.0000 

0.0000 

-  6.1929 

.5513 

0 

0 

0.0000 

0  .0000 

0.0000 

0.0000 

4# 

3233 

.9650 

0 

0 

0.0000 

0  .3000 

0.0000 

0.0000 

2. 9096 

1.3968 

96 

0 

o.oaoo 

0 .0000 

0.0000 

0.0000 

0. 

0000 

O.POOO 

97 

0 

3.7101 

3.7101 

0.0000 

0.0000 

5. 

6222 

0.0000 

0 

0 

0.0000 

0.0000 

0.0000 

0.0000 

5. 

- 5- 

8222 
*149 - 

0.0000 

.2966 

0 

0 

0 

0 

o.cooo 

0.0000 

0 .0000 
0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

5. 

6*76 

.3795 

0 

0 

0.0000 

0  .0000 

0.0000 

0.0000 

C 

-  —  -  • 

9599 

.9722 

0 

0 

0.0000 

0.0000 

0.0900 

0.0000 

4# 

0209 

.8578 

0 

0 

o.cooo 

0  .0000 

0.0000 

0.0000 

2.2619 

1.2939 

93 

0 

0.0000 

0.0000 

0.0000 

0.0000 

6  • 

0000 

0.0000 

99 

75 

-2.2899 

-2.2895 

0.0000 

75.0000 

_ ft ~nr\r\n  .  . 

0.0000 

51 

0 

3.2979 

3.2979 

0.0000  . 

(UOOOQ 

6  • 
i 

1 

2 

< 

_2. 

3 

» 

0000 

1 

3 

5 

6 

6 

.  11  . 

0.0000 

1  9 

1  9 

3  5 

3  9 

1  .  10  - 
9  19 

1  19 

53 

5 

? 

5 

1 

3 

K 

3 

75 

-2.0359 

-2  .0359 

0.0000 

75.0000 

1 

1 

1  5 

2 

2 

13 

12 

1 

1 

0 

0 

0  0 

2 

1 

2  5 

3 

12 

13 

9 

1 

1 

0 

0 

0  .  0 

3 

1 

3  9 

« 

8 

15 

26 

1 

1 

0 

0 

0  0 

_ 9_ 

_ 1  — 

9  9 

5 

26 

15 

9 

1 

1 

0. 

D  . 

-  0 _ 0. 

5 

5 

1  9 

2 

17 

29 

13 

1 

0 

0 

0  0 

b 

9 

3  S 

9 

h 

9 

26 

15 

2 

0 

0 

0  0 

7 

9 

9  5 

5 

26 

9 

* 

15 

2 

0 

0 

0  c 

6 

3 

2  9 

3 

29 

17 

6 

13 

2 

0 

0 

0  0 

9 

9 

1  10 

2 

2 

7 

2* 

17 

2 

0 

0 

0  c 

- U! 

— a — 

2  10 

3 

29 

7 

6 

17 

2 

0 

0 

0  -  0 

11 

5 

3  9 

9 

6 

17 

90 

9 

3 

0 

0 

0  0 

12 

« 

9  9 

5 

93 

IT 

9 

9 

3 

0 

0 

0  0 

13 

9 

3  10 

9 

6 

11 

90 

17 

3 

0 

0 

0  0 

19 

9 

9  10 

5 

*P 

11 

9 

17 

3 

0 

0 

0  0 

13 

10 

1  19 

2 

? 

3 

22 

7 

3 

0 

0 

0  0 

- 16 

.  12 

2  19 

3 

22 

3 

20 

7 

3 

.  0. 

0 

a _ o 

17 

10 

3  19 

9 

20 

3 

32 

11 

3 

0 

0 

0  0 

1£ 

10 

9  19 

5 

32 

3 

9 

11 

3 

0 

0 

0  0 

-1 

0. 

C 

oooo 

0  0 

0 

0 

0 

0 

0 

0 

0 

0 

0  0 

1 

1 

1  5 

2 

2 

13 

12 

1 

1 

0 

0 

0  0 

- i- 

- 1 

2..  5 

3 

12 

13 

8 

1 

1 

0 

0 

0  -0. 

3 

1 

3  9 

9 

8 

15 

26 

1 

1 

0 

0 

0  0 

9 

1 

9  9 

5 

26 

15 

* 

1 

1 

0 

0 

0  .  0 

r 

5 

1  • 

•> 

2 

17 

29 

13 

1 

0 

0 

0  0 

6 

9 

3  5 

* 

6 

9 

26 

15 

2 

0 

0 

0  0 

7 

9 

9  5 

5 

26 

9 

* 

If 

2 

0 

0 

0  0 

- 

.  .3 

2  9 

j 

29 

17 

6 

13 

2 

0 

0 

0  0 

9 

9 

1  10 
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